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ABSTRACT 

Aims. We present a spectroscopic analysis of a sample of evolved stars in M67 (turn-off, subgiant and giant stars) in order to bring 
observational constraints to evolutionary models taking into account non-standard transport processes. 

Methods. We determined the stellar parameters (Teff , log^, [Fe/H]), microturbulent and rotational velocities and. Lithium abun- 
dances (Au ) for 27 evolved stars of M67 with the spectral synthesis method based on MARCS model atmospheres. We also computed 
non-standard stellar evolution models, taking into account atomic diffusion and rotation-induced transport of angular momentum and 
chemicals that were compared with this set of homogeneous data. 

Results. The lithium abundances that we derive for the 27 stars in our sample follow a clear evolutionary pattern ranging from the 
turn-off to the Red Giant Branch. Our abundance determination confirms the well known decrease of lithium content for evolved 
stars. For the first time, we provide a consistent interpretation of both the surface rotation velocity and of the lithium abundance 
patterns observed in an homogeneous sample of TO and evolved stars of M67. We show that the lithium evolution is determined by 
the evolution of the angular momentum through rotation-induced mixing in low-mass stars, in particular for those with initial masses 
larger than 1.30 Mq when at solar metallicity. 
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1. Introduction 

The investigation of the lithium (Li) abundance in different 
stellar sites - such as open clusters - is a key element in the 
study of the chemical evolution of the Galaxy (D'Antona & 
Matteucci 1991). Regarding the formation of light elements, 
unlike beryllium (Be) and boron (B), which are formed only 
via spallation reactions involving protons or alpha particles and 
, atoms of carbon, nitrogen, and oxygen (Reeves, Fowler & Hoyle 
' 1970), lithium is mostly produced during Big Bang nucleosyn- 
I thesis. Because they are easily destroyed at low temperature 
in deep stellar layers, these light elements and especially Li 
provide strong constraints to test transport mechanisms in stellar 
interiors. 

According to standard stellar evolution, where no mech- 
anisms for the transport of matter is included in the radiative 
regions, the only episode modifying the surface abundance 
pattern of low-mass stars during their ascent of the Red Giant 
Branch (RGB) is the first dredge-up (1" DUP) occuring at the 
base of the RGB. At this stage the convective envelope deepens 
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* Based on observations collected at ESO, Paranal, Chile 
(VLT/FLAMES program ID 072.D-0309 and 074.D-0179). 
** Table 6 is only available in electronic form at 
|http : //www ■ edpsciences . org , 



and the ashes of hydrogen burning, mainly CN-cycle products, 
are dredged-up to the stellar surface, causing the decrease of 
the '^C/'^C and '^C/''*N surface ratios. The regions reached by 
the convective envelope are also lithium free and the surface 
lithium abundance therefore drops significantly. The importance 
of these variations depends on metallicity and initial mass. 

However, the large number of abundance anomalies detected 
in main-sequence and giant stars are evidence pointing to one or 
more processes of extra-mixing occurring in the stellar interiors 
(Pinsonneault, 1997). Several extra-mixing processes were sug- 
gested with the intention of explaining the Li depletion in F and 
G type stars: microscopic diffusion (Michaud 1986; Chaboyeret 
al. 1995), meridional circulation and hydrodynamical turbulent 
instabilities (Schatzman & Baglin 1991; Pinsonneault et al. 
1991; Deliyannis & Pinsonneault 1997 ; Talon & Charbonnel 
1998; Palacios et al. 2003), internal gravity waves (Garcia 
Lopez & Spruit 1991; Montalban & Schatzmann 2000; Talon 
& Charbonnel 2003). However, the physical parameters that 
control this extra-mixing are subjects of debate (Charbonnel & 
Balachandran 2000). Among the candidates, rotation-induced 
mixing alone (Palacios et al. 2003) or associated to internal 
gravity waves (Talon & Charbonnel 2003) appears to explain the 
hot and the cold side of the Li-dip observed for main-sequence 
stars in open clusters. Going beyond the turn-off, Palacios et al. 
(120031) showed that the same models that reproduce the hot side 
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of the Li-dip lead to the observed dispersion of Li abundances in 
subgiants of open clusters and of the galactic field. Concerning 
the RGB, abundance variations observed at the surface of stars 
beyond the 1" DUP indicate the existence of extra mixing 
processes at this stage (Recio-Blanco & de Laverny, 2007). 
However, rotation-induced mixing does not seem to play a major 
role in more advanced evolutionary stages, when the rotation 
decreases because of the increase of the stellar radius, as shown 
in Palacios et al. (2006) for globular cluster RGB stars. Finally 
let us mention that the second and third dredge-up episodes that 
occur during the Asymptotic Giant Branch phase once again 
modify the surface abundance patterns, and that in these phases 
abundance anomalies also indicate the action of an extra-mixing 
process that is not yet clearly identified. Li depletion also seems 
to be related to the stellar age (Herbig 1965; Duncan 1981; 
Soderblom 1983; Fekel and Balachandran 1993), and it is also 
sometimes considered to be connected with other parameters 
such as metallicity, activity, and mass (Pallavicini et al. 1987; 
Spite & Spite 1982; Randich et al. 1994). If the Li depletion is 
related to the stellar age, one could expect a correlation between 
lithium abundance and rotational velocity for stars of the same 
mass, same metallicity and same spectral type (Skumanich 
1972). In agreement with this, Zahn (1992) and Pinsonneault et 
al. (1990) postulated that the Li depletion in late-type stars is 
directly related to the loss of angular momentum. 

The age of clusters and consenquently that of their individ- 
ual stars can be determined from isochrone fitting. Cluster stars 
may therefore more efficiently serve for the study of changes in 
mixing-sensitive abundances and for the search of anomalies in 
chemical compositions. Open clusters are moreover important 
laboratories, from the study of which one can address the 
question of Li evolution in stellar interiors and surfaces, because 
they contain a significant sample of stars in a wide range of 
mass with the same origin and same initial chemical compo- 
sition. In the last decades several observational diagnostics 
have discovered element abundances of stars belonging to 
open clusters of different ages with the help of high-resolution 
spectroscopy, and even more recently with the multi-objects 
instruments. The evolution of Li and Be abundances along 
specific parts of the colour-magnitude-diagrams (CMD) of open 
clusters (main-sequence, turn-off (hereafter TO), subgiant and 
red giant branch) has been investigated to test the physics and 
the extra-mixing processes included in stellar evolution models. 
Randich et al. (2002 & 2007) have investigated both Be and Li 
abundances in late-F and early-G main-sequence and subgiant 
stars from five different open clusters with ages ranging from 
50 Myr up to 4.5 Gyr. They confirmed shallow mixing in 
stellar interiors, which is able to transport surface material deep 
enough for Li burning to occur, but not deep enough for Be 
burning. Very recently, Smiljanic et al. (2010) have analysed 
Be and Li abundances along the main sequence and RGB in 
the open cluster IC 4651 (with an age estimated around 1.5 
Gyr). They show that both the evolution and dispersion of Be 
and Li abundances are successfully reproduced by theoretical 
predictions from hydrodynamical models taking into account 
non-standard mixing processes (rotation-induced mixing, 
atomic diffusion, and thermohaline mixing). 

The present study focusses on evolved stars of the open 
cluster M67 (- NGC 2682), which has been used for the 
past three decades as an important laboratory for studying 
stellar evolution (Burstein et al. 1986; Carraro et al. 1996, and 
references therein). M67 is commonly mentioned as a solar 



age cluster, i.e. one of the oldest open clusters of the Galaxy. 
However, its age is still a matter of debate, and estimates vary 
from 3.5 to 4.5 Gyr A recent estimate (Sarajedini et al., 2009) 
based on the comparison of a M67 proper-motion-cleaned 
colour-magnitude diagram to theoretical isochrones points to 
low values around 3.5 - 4.0 Gyr, confirming previous results 
from VandenBerg & Stetson (2004) or Michaud et al. (2004), 
who used evolutionary models taking into account microscopic 
processes such as atomic diffusion. On the other hand, spec- 
troscopic studies devoted to stars from M67 usually refer to a 
higher age of about 4.5 Gyr (e.g., Randich et al. 2006). Using 
spectroscopy and photometry, several authors determined an 
average iron abundance [Fe/H] for M67 very close to the solar 
value. 

Previous Li observations in M67 stars show that a real 
dispersion in the lithium abundance of main-sequence objects 
exists even though at this phase the convective envelope is 
too superficial to reach the layer where Li destruction can 
occur (Pasquini et al. 1997; Jones et al. 1999; Randich et al. 
2002; 2007). While standard models appeared inefficient to 
explain this Li dispersion observed in unevolved M67 stars, the 
rotational mixing seems to be a good candidate (Pasquini et al. 
1997). Considering M67 subgiants. Sills & DeHyannis (2000) 
also find that stellar models with slow mixing caused by rotation 
are the most suitable to explain the evolution of low-mass 
stars. By measuring rotational velocities for 28 stars from the 
main-sequence to the giant branch in M67, Melo et al. (2001) 
provided an analysis of the history of the angular momentum 
of stars of 1.2Mq. It verifies that these velocities probably obey 
different laws for the evolution of the angular momentum on the 
main-sequence and on the giant branch. 

Our observations of M67 stars along an evolutionary se- 
quence from the TO to the top of the RGB offers the oppor- 
tunity to shed light on the interplay between Li and rotation 
evolution and to test the extra-mixing scenario in stellar struc- 
tures. In Section 2 we present the observational data and our 
method to determine stellar parameters (effective temperature, 
surface gravity, metallicity), Li abundances, and rotational ve- 
locities (v sin /). In Section 3 we present the stellar models we 
used to study the Li evolution in evolved stars of M67, and our 
conclusions are drawn in Section 4. 



2. Spectroscopic observations and data analysis 

Our study is based on a sample of 28 post-main-sequence stars 
of the open cluster M67. Figure [T] displays a colour-magnitude 
diagram of the open cluster M67 (photometry from Montgomery 
et al., 1993). We overplotted a best-fitting isochrone of 3.7 Gyr 
provided by VandenBerg (private communication) based on the 
models wiffi diffusion by Michaud et al. (2004). 

These stars are divided into three groups regarding their evo- 
lutionary stages: (i) TO stars, (ii) subgiant branch stars, and (iii) 
RGB stars, including three clump giants. According to a stan- 
dard evolutionary scenario, TO and subgiant stars must present 
a fast expansion of their convective envelopes, while giant stars 
must have passed akeady through the first dredge-up episode, 
bringing to the surface CNO products, hence modifying the sur- 
face abundances. Let us mention that among our sample is the 
Li -rich subgiant S1242, which we analysed in detail in Canto 
Martins et al. (2Q06j). 
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Fig. 1. Colour-magnitude diagram of the open cluster M67. (V, 
B - V) photometry is taken from Montgomery et al. (1993) and 
a colour excess of E{B - V) - 0.05 has been assumed, as in 
Melo et al. (2001). Stars from our sample are shown as open dot 
circles. The track is an isochrone at 3.7 Gyr from Vandenberg 
(private communication) based on models with diffusion from 
Michaud et al. (2004). 



2.1. Spectroscopic observations 

The observations of these stars were carried out with the 
VLT/Unit2 ESO telescope (Paranal, Chile), using the FLAMES- 
UVES spectrograph (Pasquini et al. 2002), during January 
2004, January 2005 and March 2005. Using photometry and 
astrometry data of M67 from EIS pre-FLAMES (Momany et al. 
2001) we selected our stars near the centre field with a = 8 : 48 
and (5 = +11 : 46. We used two different fiber configurations. 
For the brightest stars (V < 12, i.e., RGB stars) we observed 
four stars, with four fibers dedicated to the sky with the exposure 
time set to 1 500 sec, providing a mean signal-to-noise ratio 
(S/N) of 50 (per pixel), over the lithium region (6 700 - 6 720 
A). For the faintest stars (V > 12, TO and subgiant stars), seven 
fibers were used for the stars and one for the sky. The exposure 
time for these groups was set to 1 500, 2 580, and 3 000 sec, 
providing a S/N between 65 and 100 (per pixel). The stellar 
sample was observed with the UVES red arm centred at 580 nm, 
covering a wavelength range of 420-680 nm, and also centred 
at 860 nm, covering a larger wavelength range of 660-1060 
nm. The observations have a resolving power of R ~ 47000 (1 
arcsec sky aperture). Moreover, some stars (marked in Table [T] 
with a star) we only observed with the UVES red arm centred at 
580 nm. 

The spectra were reduced with the standard FLAMES/UVES 
data reduction pipeline following the usual steps of reduction 
(bias, flat-field, and background corrections, order fiber defini- 
tion, wavelength calibration of the spectra with a ThAr spectrum 
and extraction of the spectra). Concerning the sky subtraction, 
for bright stars, four fibers were allocated for the sky on each 



plate. The average of the sky fibers was subtracted from each 
star spectrum. However, for the faintest stars only one fiber was 
dedicated to the sky and then subtracted from the star spectrum. 
Then we used IRAF to normalize them to a pseudo continuum 
and to bring the reduced spectra to rest. The radial velocities 
for our stellar sample were calculated with the IRAF fxcor 
task. The stellar spectra were cross-correlated with a spectrum 
of the Sun (Hinkle et al. 2000). The shifts were then computed 
into radial velocities of the stars and heliocentric correction was 
applied. The results are presented in Table 1 and the errors in 
these velocities show a stability at a level of a few hundreds of 
na/s. 

The probability that these stars belong to this stellar cluster 
is larger than 70%, according to Sanders (1977) and their esti- 
mated radial velocities should be consistent with the estimated 
radial velocity for M67 of 33.0 km s^' (Friel & Janes 1993). 
This is not the case for the star SI 000 with 0% of membership 
probability and a radial velocity of 42.8 km s"', which is very 
different from the group of stars and is therefore rejected. 

Finally, we also combined different exposures collected at 
different dates for some TO and subgiant stars, in order to reach 
a good quality of the data and to increase the mean S/N over the 
Li region to 100 per pixel. 

A complete log of observations is presented in Table [1] Star 
identifiers are from Sanders (1977). The radial velocity correc- 
tion that was applied is mentioned in the last column (which 
agrees very well with the values produced by Melo et al. 2001). 

2.2. Spectral syntiiesis 

We performed spectral synthesis analysis on our 27 sample 
stars, to derive stellar parameters (Tgff , logg , ^, [Fe/H]) and to 
measure rotational velocities and lithium abundances. We used 
the MARCS models of stellar atmospheres (Gustafsson et al. 
2008), which are based on plan-parallel and spherical models at 
local thermodynamical equilibrium (LTE). The turbospectrum 
spectral synthesis tools (Alvarez & Plez 1998) and interpolation 
routines on model atmospheres were used. Solar abundances 
were taken from Asplund, Grevesse & Sauval (2005) and the 
collisional damping treatment was performed based on the work 
of Barklem et al. (Barklem et al. 2000a,b; Barklem & Piskunov 
2003, Barklem & Aspelund-Johansson 2005). To compute 
synthetic spectra, we took taken into account atomic (see below) 
and molecular line lists: TiO (Plez 1998), VO (Alvarez & Plez, 
1998), and CN and CH (Hill et al. 2002). 

In order to improve the precision on the atmospheric 
parameters and on the Li abundances (Au - tog(^) + 12) 
presented in Canto Martins et al. (2007), we analysed our 
sample using the curve-of-growth method, which is based on the 
measurement of equivalent widths (EWs) of a large number of 
iron lines present in our FLAMES-UVES observations. Indeed, 
weak features (mainly caused by the Fe I line at 6707.44 
A) are known to blend the Li line (at 6 707.78 A), which 
is used to derive Au ■ Hence, a very good precision on the 
estimate of the Fe abundance and on all atmospheric parameters 
is required to ensure a determination of Au with a high accuracy. 

To do so, we first calibrated the oscillator strength values 
(logg/) of 92 Fe I and 14 Fe II spectral lines issued from the 
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Table 1. Log of the spectroscopic observations of our 28 sample stars. The stars were observed using two different setups, except 
for those marked with (★), which were only observed using the setup centred at 580 nm (see text). 



Sanders 


V 


(B-V) 


0^2000 


1^2000 


Julian Day 


Exposure 


Vrad 


ID 










+ 2453000 


time (sec) 


(km.s"') 


S774 


12.93 


0.85 


08:50:49.93 


11:49:12.90 


009.74 


3000 


33.5 












011.77 


1500 


34.2 












031.76 


3000 


34.0 












032.76 


3000 


33.6 












033.64 


3000 


34.0 


S806 


12.78 


0.81 


08:51:00.17 


11:54:32.10 


009.74 


3000 


33.4 












011.77 


1500 


34.3 












031.76 


3000 


34.0 












032.76 


3000 


33.6 












033.64 


3000 


34.0 


S978 


9.72 


1.37 


08-51 17 48 


11:45:22.70 


008.77 


1500 


34.7 












008.79 


1500 


34.6 


SIOOO* 


12.84 


0.78 


08-5 1-23 83 


11:47:16.20 


382.73 


2580 


42.8 












382.76 


2580 


42.9 


SlOlO 


10.48 


1.11 


08:51:22.81 


11:48:01.80 


008.77 


1500 


34.1 












008.79 


1500 


33.8 


S1016 


10.3 


1.26 


08:51:17.11 


11:48:16.30 


008.77 


1500 


34.6 












008.79 


1500 


34.4 


S1034 


12.65 


0.61 


08:51:18.55 


11:49:21.70 


028.74 


3000 


34.6 












028.78 


3000 


34.6 












033.69 


3000 


34.8 












033.73 


3000 


34.9 


S1074 


10.59 


1.12 


08:51:12.68 


11:52:42.20 


007.76 


1500 


34.2 












007.78 


1500 


34.7 


S1231 


12.93 


0.92 


08:51:29.38 


11:45:27.80 


009.74 


3000 


32.9 












011.77 


1500 


33.6 












031.76 


3000 


33.5 












032.76 


3000 


33.1 












033.64 


3000 


33.5 


SI 239 


12.75 


0.76 


08-51-44 03 


11:46:24.60 


028.74 


3000 


32.8 












028.78 


3000 


32.9 












033.69 


3000 


33.2 












033.73 


3000 


33.2 


S1242 


12.72 


0.68 


08:51:36.04 


11:46:33.70 


009.74 


3000 


37.5 












011.77 


1500 


38.4 












031.76 


3000 


34.1 












032.76 


3000 


34.4 












033.64 


3000 


35.1 


SI 245 


12.93 


0.92 


08:51:44.77 


11:46:46.20 


028.74 


3000 


33.1 












028.78 


3000 


33.1 












033.69 


3000 


33.3 












033.73 


3000 


33.4 


SI 254 


11.52 


1.05 


08:51:45.08 


11:47:46.00 


007.76 


1500 


32.9 












007.78 


1500 


33.5 


S1268 


12.65 


0.58 


08:51:49.97 


11:49:31.40 


009.74 


3000 


32.7 












011.77 


1500 


33.4 












031.76 


3000 


33.2 












032.76 


3000 


32.9 












033.64 


3000 


33.3 


SI 273 


12.22 


0.57 


08-51-39 26 


11-50 04 00 


028.74 


3000 


34.6 












028.78 


3000 


34.6 












033.69 


3000 


34.7 












033.73 


3000 


34.8 


SI 275 


12.56 


0.59 


08:51:37.43 


11:50:05.40 


009.74 


3000 


33.4 












011.77 


1500 


34.1 












031.76 


3000 


33.9 












032.76 


3000 


33.7 












033.64 


3000 


34.0 


SI 277 


11.63 


1.05 


08:51:42.35 


11:50:07.80 


008.77 


1500 


34.6 












008.79 


1500 


34.3 


S1279 


10.55 


1.12 


08:51:28.98 


11:50:33.00 


007.76 


1500 


33.3 












007.78 


1500 


33.9 


S1288 


11.33 


1.07 


08:51:42.37 


11:51:23.10 


007.76 


1500 


33.4 












007.78 


1500 


34.0 


S1293 


12.09 


1.01 


08:51:39.42 


11:51:45.90 


028.74 


3000 


34.2 












028.78 


3000 


34.2 












033.69 


3000 


34.3 












033.73 


3000 


34.5 
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Sanders 


V 


(B-V) 


Q'2000 


(52000 


Julian Day 


Exposure 


I'rad 


ID 










+ 2453000 


time (sec) 


(km.s"' ) 


S1305* 


12.23 


0.99 


08:51:35.79 


11:53:35.00 


398.79 


2580 


34.3 












431.68 


2580 


34.2 


S1319* 


12.91 


0.97 


/^O C i Act Ctf\ 

08:51:48.80 


11:56:52.0 


398.79 


2580 


34.4 












431.68 


2580 


34.3 


SI 323 


12.83 


0.78 


08:51:35.43 


11:57:56.80 


028.74 


3000 


33.2 












028.78 


3000 


33.2 












033.69 


3000 


33.5 












033.73 


3000 


33.5 


S1438* 


12.89 


0.82 


08:52:11.35 


11 AC C f\ 

11:45:37.50 


398.79 


2580 


33.1 












431.68 


2580 


33.0 


S1487 


12.63 


0.64 


f\Ct C f\ A O 1 

08:52:04.81 


11:58:28.80 


398.79 


2580 


33.4 












431.68 


2580 


33.3 


S1607* 


12.62 


0.56 


08:52:21.45 


11:50:41.70 


398.79 


2580 


33.7 












431.68 


2580 


33.7 


S2207* 


12.65 


0.62 


08:51:32.47 


11:47:52.50 


398.79 


2580 


31.9 












431.68 


2580 


31.9 


S2208 


12.78 


0.82 


08:51:32.42 


11:48:01.30 


028.74 


3000 


33.0 












028.78 


3000 


33.0 












033.69 


3000 


33.1 












033.73 


3000 


33.1 



VALD database (Kupka et al. 1999). These lines appear reason- 
ably unblended throughout the spectral domain (4 200-8 000 A) 
of the high-resolution spectra of the Sun and of Arcturus 
(Hinkle et al. 2000). For each line the central wavelength, 
excitation potential (Xexc), and our value on log gf improved by 
an inverse solar and Arcturus analysis are presented in Table |6] 
(online material). The equivalent widths (EW) from the Sun and 
from Arcturus spectra (Hinkle et al. 2000) are also reported. 
They were measured with an automatic procedure using the 
DAOSPEC package (Stetson & Pancino 2008). In order to 
check our corrections on log gf values, we used our measured 
EWs to derive the Fel and Fell abundances for the Sun and for 
Arcturus. For the Sun, we adopted the following atmospheric 
parameters: Teff= 5777 K (Neckel 1986), ^ = 1.0 km s"' 
(Ruedi et al. 1997), \ogg = 4.44 (Allen 1973); and for Arcturus: 
Teft =4300 K, ^ = 1.6 km s-i , \ogg = 1.8 (ZoccaU et al. 2004). 
The solar iron abundances we derived, log n(Fel) -1A9 + 0.03 
dex and logn{FeII) = 7.48 + 0.04 dex, agree well with the 
values of Asplund, Grevesse & Sauval (2005) and also with 
the value of the Fe II abundance found by Biemont et al. 
(1991). For Arcturus we found log n{FeI) - 7.02 + 0.05 dex and 
log n{FeII) = 7.05+0.06 dex, which agrees with determinations 
from Peterson et al. (1993) and from Carraro et al. (2004). 



lines for the Li region. 



: (a) 

1^ -J 






. 



o 7.4 
7.3 



(b) 




Then we estimated for our 27 sample stars the effective 
temperatures, microturbulent velocities, and stellar surface 
gravities by imposing an excitation equilibrium for the Fe I 
abundances, an equilibrium between the Fe I abundances and 
the equivalent widths, and a Fel/Fell ionization equilibrium as 
presented in Fig. |2]for the TO star S1275. The initial values 
adopted for the atmospheric parameters were the ones published 
in Canto Martins et al. (2007). In this previous work, the 
TefF were estimated from photometry (Montgomery et al. 1993; 
Houdashelt et al. 2000), from the H„ and lines and also the 
Fe I lines ratio in the Li region. The log g was estimated using 
evolutionary models. For the metallicities, the solar metallicity 
was adopted as a first guess and then adjusted by fitting the Fe I 



log (EWA) 

Fig. 2. Analysis of the turn-off star S 1275 : (a) Fe lines excitation 
equilibrium, (b) relations between Fe abundances and equivalent 
widths. Open and filled circles represent the Fe I and Fe II lines, 
respectively. The solid line represents the mean iron abundance 
and the dashed lines represent the Icr of the distribution. 



With this method the metallicities [Fe/H] were calculated 
for a solar abundance log n(Fe)Q = 7.49 (as derived here) 
and the [Fe/H] errors are based on the standard deviations in 
the Fe I lines abundances only. The accuracies obtained on 
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atmospheric parameters are +70 K for Tefi , ±0.2 dex for log g 
and +0.2 km s"' for ^. 

Table |2] provides the atmospheric parameters adopted to 
compute the synthetic spectra used to derive the projected 
rotational velocity (v sin /) and Au for each star in our sample. 
To derive the projected rotational velocity, we used the same 
procedure as in de Medeiros et al. (2006). We convolved the 
resulting spectra (taking into account the instrumental profile 
of FLAMES-UVES) with a rotational profile to adjust the 
broadening observed on a group of Fel lines in the Li region 
(6700-6720 A). The vsini values we measured (presented in 
Table|2| have an accuracy of +1.0 km s"'. 

Finally, the determination of Au was made by the best fit 
between observed and synthetic spectra, obtained on the Li line 
at 6707.78 A. For each star, the total error on the Li abundance 
was estimated by computing the quadratic sum of errors induced 
by errors on individual parameters (Teft , log g , ^, [Fe/H], v sin /) 
and also errors associated with the best-fit determinations. As 
discussed by De Laverny et al. (2003) and Lebre et al. (2006), 
the main source of uncertainty in the abundances determina- 
tions is caused by errors in the Tefj determination. They have 
estimated that Tefj measurements with errors smaller than 200 K 
lead to errors smaller than 0.2 dex in the derived metallicity 
and lithium abundance. We made some tests, and the obtained 
errors are about the same as proposed in these previous studies. 
However, we found that errors in the parameters \ogg and ^ 
also influence the determination of the metallicity of the stars in 
our sample. We found that the errors obtained in [Fe/H] caused 
by errors in \ogg are smaller than those measured when we 
analyse the errors in the metallicity caused by ^. Taking into 
account the two sources of errors, we obtained a maximum error 
of 0.16 dex in [Fe/H]. In its last column. Table |2] also presents 
ALi values and their accuracies. As an illustration we show in 
Fig. |3]a zoom on the Li region together with synthetic spectra 
for several values of Au for the TO star S1275. 

The stellar parameters, microturbulent and rotational veloc- 
ities and Au measured in this study are more accurate than the 
values produced by Canto Martins et al. (2007) because of the 
large number of Fe I and Fe II lines used in the characterization 
of the stellar sample. 

2.3. Comparison with otiier spectroscopic studies 

We compared the parameters we determined with previous 
spectroscopic studies available in the literature. Concerning the 
metallicity of M67 we found from our sample stars a mean 
metallicity for our evolved stars of [Fe/H] - -0.05 + 0.04. This 
value agrees very well with the lowest values from the literature 
as [Fe/H] = -0.05 (Canterna et al. 1986), [Fe/H] = -0.07 
(Anthony-Twarog 1987), [Fe/H] = -0.08 (Friel & Jane 
1991), [Fe/H] = -0.09 (Friel & Jane 1993), [Fe/H] = -0.05 
(Shetrone & Sandquist 2000) and [Fe/H] = -0.03 (Tautvaisiene 
et al. 2000). However, other studies point to higher values for the 
mean metalHcity of M67 as [Fe/H] - 0.06 (Nissen et al. 1987), 
[Fe/H] = 0.04 (Garcias Lopez et al. 1988), [Fe/H] = 0.02 
(Friel & Boesgaard 1992), and recently, [Fe/H] = 0.02 + 0.14 
(Yong et al. 2005), [Fe/H] = 0.03 + 0.01 (Randich et al. 2006, 
hereafter R06) and [Fe/H] = 0.05 + 0.02 (Pancino et al. 2010). 



Table 2. Estimated atmospheric parameters for our sample stars. 
Rotational velocities and lithium abundances derived from our 
spectral synthesis analysis are also given. The identifiers (ID) 
are from Sanders (1977). 



ID 


Teff 






[Fe/H] 


V sin i 






(K) 


(dex)(km.s 


) (dex) (km.s^' 


) 


Turn-off stars 


SI 273 


6159 


4.01 


1.86 


-0.04±0.04 


8.0 


1.90+0.11 


SI 607 


6127 


3.81 


1.77 


-0.11+0.06 


3.5 


1.70+0.10 


SI 275 


6050 


4.00 


1.64 


-0.04+0.04 


4.9 


2.15+0.11 


S1034 


6020 


3.90 


1.94 


-0.08+0.03 


4.0 


1.30+0.11 


S2207 


6000 


3.90 


1.62 


-0.03+0.03 


4.0 


1.20+0.12 


S1268 


5996 


3.86 


1.75 


-0.11+0.07 


3.5 


0.90+0.12 


S1487 


5940 


3.81 


1.67 


-0.05+0.03 


2.5 


1.10+0.13 


Subgiant stars 


SI 242 


5810 


3.90 


1.50 


-0.04+0.06 


6.1 


2.70+0.13 


S1323 


5654 


3.90 


1.50 


+0.03+0.04 


3.3 


<1.00 


S1239 


5644 


3.80 


1.41 


+0.00+0.03 


3.1 


<0.80 


S806 


5461 


3.80 


1.25 


+0.08+0.03 


4.1 


0.00+0.19 


S2208 


5429 


3.90 


1.24 


+0.01+0.03 


3.6 


0.70+0.14 


S1438 


5420 


3.80 


1.27 


-0.06+0.03 


2.2 


<0.00 


S774 


5240 


3.70 


1.20 


-0.08+0.03 


3.2 


<0.00 


Giant stars 


SI 245 


5170 


3.61 


1.19 


-0.02+0.03 


3.0 


0.00+0.15 


S1231 


5130 


3.60 


1.12 


-0.02+0.03 


3.0 


<0.10 


S1319 


5104 


3.61 


1.16 


-0.07+0.03 


3.0 


<0.40 


S1293 


4970 


3.30 


1.32 


-0.01+0.03 


3.2 


-0.20+0.22 


S1305 


4940 


3.20 


1.18 


-0.08+0.03 


2.8 


< -0.20 


SI 254 


4820 


2.91 


1.30 


-0.03+0.04 


3.0 


< -0.40 


S1277 


4820 


3.00 


1.26 


+0.01+0.05 


3.0 


< -0.40 


SI 279 


4779 


2.72 


1.57 


-0.01+0.08 


2.5 


-0.50+0.23 


S1288 


4773 


2.90 


1.32 


-0.01+0.04 


2.5 


-0.20+0.21 


S1074 


4750 


2.60 


1.62 


-0.07+0.04 


2.5 


-0.40+0.24 


SlOlO 


4748 


2.60 


1.58 


-0.03+0.07 


2.5 


<0.00 


S1016 


4430 


2.31 


1.52 


-0.05+0.05 


2.0 


-0.50+0.25 


S978 


4260 


1.90 


1.63 


-0.15+0.06 


2.0 


-1.00+0.27 



Table 3. Comparison for the star SlOlO between atmospheric 
parameters from Tautvaisiene et al. (2000) (TOO), Young et al. 
(2006) (Y05), and Pancino et al. (2010) (PIO) and those derived 
in the present work. 



Parameters 


TOO 


¥05 


PIO 


Our work 


T.ff (K) 


4730+100 


4700+100 


4650+100 


4748+70 




2.4+0.3 


2.3+0.3 


2.8+0.2 


2.6+0.20 


^(km s-') 


1.6+0.3 


1.34+0.2 


1.3+0.1 


1.58+0.20 


[Fe/H] 


-0.01+0.11 


0.00+0.10 


0.06+0.01 


-0.03+0.07 



Comparing some stars in our sample with other studies in 
greater detail, we have the star SlOlO that was also studied 
by Tautvaisiene et al. (2000), Yong et al. (2005), and Pancino 
et al. (2010). Table [3] compares the values of the atmospheric 
parameters and metallicity measured for this star Evidently 
the atmospheric parameters obtained in all studies agree well. 
However, we measured a lower metallicity compared to their 
results. One possible explanation for this is the higher solar 
iron abundance adopted by these studies (logn{FeI) - 7.54 
for Yong et al. 2006 and log n(F el) = 7.50 for Pancino et 
al. 2010). Another explanation is that both the superficial 
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Table 4. Comparison for SI 034 and SI 239 between atmospheric to compute the model atmosphere and synthetic spectra, 
parameters from Randich et al. (2006) and lithium abundances 
from Randich et al. (2007) and those derived in the present work. 





S1034 


S1239 


Randich 


Our work 


Randich 


Our work 


Teff(K) 


5969±70 


6020+70 


5477+70 


5644+70 




4.0+0.25 


3.9+0.20 


3.8+0.25 


3.8+0.20 


^(km s-') 


1.5+0.15 


1.94+0.20 


1.25+0.15 


141+0.20 


[FelH] 


+0.01±0.05 


-0.08+0.03 


+0.02+0.03 


+0.00+0.03 


Au 


<1.14 


1.30+0.11 


1.12+0.10 


<0.80 




6704 6706 6708 6710 

Wavelength { A) 

Fig. 3. Lithium line spectroscopic region of the turn-off star 
SI 275: (a) observation (dotted line) and synthetic spectra com- 
puted with the atmospheric parameters listed in Table|2] and with 
Au = 2.10, 2.15, 2.20 dex (dashed, solid and dashed-dot lines re- 
spectively). The best synthetic spectrum is clearly the one com- 
puted with Au - 2.15. (b) Residual difference between synthetic 
and observed spectra. 



gravity and the microturbulent velocity affect the metallicity 
of the star. In this case, as the value of f obtained in our 
study is higher than that obtained by Yong et al. (2005) 
and Pancini et al. (2010), this difference leads to errors on 
[Fe/H] smaller than 0.12 dex in SIOIO. This difference may 
also explain the low mean metallicity of M67 found in our study. 

Two other targets, S1034 and S1239 (Fig. |4|l, were also 
included in the spectroscopic study of Randich et al. (2006). For 
these two objects. Table |4]provides the comparison of the stellar 
parameters between R06 and the present study. Both analyses 
appear to agree fairly well within the error bars for stellar 
parameters determination, even if we derive a lower metallicity 
for S1034 (but still consistent within the error bars) and a 
higher microturbulent velocity. However, the Au values present 
differences that may be explained by the adopted methodology 



S1034 



y — 



S1239 - 



6706 6708 
Wavelength ( A) 



Fig. 4. Li region of the stars SI 034 and S1239, which were also 
studied by Randich et al. (2006). Observations are represented 
with a dotted line and synthetic spectra by a dashed line. The 
vertical line indicates the Li feature at 6707.78 A. 



Finally, our rotational velocities (fsin/) agree well with 
the results obtained by Melo et al. (2001) on evolved stars of 
M67. Moreover our study statistically extends the number of 
M67 evolved stars with measured v sin /. In this way, our efforts 
complement the existing values of atmospheric parameters and 
V sin / values for M67 objects, and can assist studies that are 
mainly dedicated to the evolution of the angular momentum in 
the open cluster M67. 

Up to date, as mentioned in different studies (Pilachowski et 
al. 1988, Balachandran 1995), subgiant and giant stars in M67 
show no significant Ali , which points to a severe Li depletion 
after the TO. For unevolved stars of M67, Jones et al. (1999) 
and Randich et al. (2002, 2007) show a large dispersion on the 
Ali . Balachandran (1995) also studied the Li abundance in M67 
evolved stars, but only upper limits were estimated. In this con- 
text, the analysis we performed on our sample complements the 
previous observational studies mainly devoted to stars of M67 
near the TO. Our data present the same dispersion as established 
by Jones et al. (1999) and Randich et al. (2002 and 2007) and 
shows a gradual decrease of Ali for stars cooler than 5500 K 
(see Fig.|71 next section). The present result points to a depletion 
of Li in excess to standard predictions as pointed out by Pasquini 
et ak (1997). 



3. Connecting the Li abundances and the rotational 
velocities 

The lithium abundances derived for our stars clearly indicate that 
this element has been depleted at their surface already during the 
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main sequence evolution. In an attempt to consistently explain 
the evolution of both the lithium abundances and the rotational 
velocities, we computed a grid of stellar evolution models with 
rotation and compared the predictions to our sample. 



more detailed treatment for microscopic diffusion and radiative 
forces. Bellow we will therefore consider an age of 3.7 Gyr for 
M67 and all discussed predictions will refer to values obtained 
when the models reach this age. 



3.1. Stellar evolution models with rotation 



We computed a grid of stellar evolution models in the mass 
range 1.30 Mq to 1.37 Mq (see Table |5]l including rotation, 
microscopic diffusion, and thermohaline mixing (as described 
in Charbonnel & Zahn 2007) using STAREVOL V3.00 (Siess 
12006 : Siess et al. 2008; Decressin et al. 2009 ). We used the so- 
lar abundances given by Asplund et al. (2005) as a reference 
for our stellar evolution models, and the corresponding opacities 
were generated with the OPAL Opacity code (Iglesias & Rogers 
[1996Q. Convective zones are defined using the Schwarzschild 
criterion and the mixing length formalism with an adopted mix- 
ing length parameter aMLT = 1-78 that was calibrated on the 
Sun. The models have a metallicity of [Fe/H] = -0.03, which 
translates into a metal content Z = 0.01148. The initial helium 
mass fraction Yini = 0.2688 was calibrated for the solar model. 
The initial lithium content is set to Au = 3.2CB Mass loss is 
included in our models from the ZAMS and beyond following 
Reimers (119751 ). Microscopic diffusion (gravitational settling) of 
chemical species is accounted for using the approximation of 
Paquette et al. (1986) for the microscopic diffusion coefficient 
and the expressions given by Montmerle & Michaud (119761 ) for 
the microscopic diffusion velocity. We adopt the Maeder & Zahn 
( 119981 ) formalism to account for the transport of angular mo- 
mentum and of chemical species by meridional circulation and 
shear-induced turbulence in the models. The formalism and its 
introduction in STAREVOL are described at length in Palacios et 
al. (120031 1200 6^. and we refer the reader to these papers for fur- 
ther details. Let us nonetheless recall that we use the prescription 
from Zahn ( 119921) for the horizontal turbulent diffusion coeffi- 
cient Di, with Ct, - 1 .0. We also assume that convective regions 
rotate as solid bodies. For those models with equatorial velocity 
at the ZAMS of about 19 km.s"', we assume angular momen- 
tum losses associated with a magnetic torque as described by 
Kawaler (IT988b : 



(Q. < Q.sat) 



[dt , 



R 



Rc 



1/2 



— = -KQ.n,„2 — — (O > 



M 



Mr 



-1/2 



(1) 



(2) 



The parameter K in Kawaler's law is related to the magnitude 
of the magnetic field strength, and Q,,,, = 14Q.Q following 
Bouvier et al. (119971 ). Lacking observational data on the ZAMS 
to constrain the evolution of angular momentum of our sample 
stars, we use the vsin/ derived by Melo et al. (2001) for main- 
sequence stars of M67 as a guideline to calibrate the braking. 
For masses 1.30 to 1.37 Mq we adopt K = 8x10^''. 

We checked that a model of 1.25 M© including microscopic 
diffiision, no rotation, and a mixing length parameter calibrated 
for the Sun reaches the turn-off at 3.75 Gyr, in fair agreement 
with the value obtained by Michaud et al. (120041 ). who used a 



' The tables were generated using the web interface available at 
lhttp://opalopacity.llnl.gov/ 

This abundance is obtained using the solar meteoritic abundance 
given by Aslpund et al. J2005I ) scaled down to the metallicity of M67. 




7000 



6000 



5000 



4000 



Fig. 5. log g versus Tefj for our sample of stars, represented by 
black circles. Values at 3.7 Gyr for the rotating models with and 
without angular momentum extraction are represented by filled 
red pentagons and filled green squares respectively. The dotted 
line indicates the turn-off (TO) position and the location of the 
first DUP is also indicated. 



Table |5] presents the results of our modelling in terms of 
effective temperature and gravity, luminosity, equatorial veloc- 
ity, and lithium abundance at the ZAMS, the turn-off (e.g. TO) 
and 3.7 Gyr (e.g. M67). Models with and without braking have 
an initial angular momentum of about 5.2x10'*'^ and 3.2x10"^^ 
g.cm^.s"' respectively. 

3.2. Evolutionary status and rotation 

Our choice of masses for the models computed are based on 
the 1.3 Mq turn-off mass suggested by Michaud et al. (120041) 
for their 3.7 Gyr best-fitting isochrone to the M67 colour- 
magnitude diagram, and taking into account that our sample 
includes subgiant and giant stars that should be more massive. 
Figure |5] shows the surface gravity as a function of temperature 
for our sample stars together with our models at the age of 3.7 
Gyr. Both types of models, whether they are slow or moderate 
rotators on the ZAMS, fall among the data points. The data are 
well reproduced by rotating models that have passed the turn-off 
(see Table|5]), and that have initial masses between 1.30 Mq and 
1.37 Mo. 

Figure |6] presents the rotation velocity v sin / of our data 
versus temperature and overplotted are the coiTesponding 
V sin / from our models at 3.7 Gyr. The models have an actual 
equatorial rotation velocity higher than the values plotted 
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Table 5. Characteristics of stellar evolution models with rotation. Effective temperature, surface gravity, luminosity, surface rotation 
velocity, lithium abundance, age, and associated evolutionary phase are given for models with initial masses between 1.30 Mq and 
1.37 Mq. In the last column, "M67" refers to the models predictions at the adopted age for the cluster of 3.7 Gyr For each quantity, 
the values for models without (/Tk - 0) and with (/Tk ^ 0) magnetic braking on the ZAMS are indicated, the latter is shown in bold 
face. 



mass Teff log g L Y,^,f Au age phase 
(Me) (K) (dex) (Lo) (tan/s) (dex) (Gyr) 





^K=0 


A-K ^0 


Kk = 


A-K ^0 


Kv. = Q 




Kk = Q 




Kk = 


Kk*0 


J^K = 


Kk*0 




1.30 


6699 


6701 


4.37 


4.37 


2.75 


2.75 


12.05 


18.9 


3.18 


3.18 


4.48 10-- 


A.16 10-^ 


ZAMS 




6439 


6440 


3.99 


3.99 


5.55 


5.57 


10.86 


8.88 


2.99 


1.26 


3.36 


3.38 


TO 




6212 


6218 


3.89 


3.89 


6.11 


6.15 


7.97 


8.03 


2.93 


1.14 


3.7 


3.7 


M67 


1.31 


6729 


6729 


4.36 


4.36 


2.85 


2.85 


11.92 


18.87 


3.18 


3.18 


4.50 10-^ 


4.50 10-^ 


ZAMS 




6457 


6472 


3.99 


3.99 


5.70 


5.70 


10.72 


9.00 


2.93 


1.27 


3.24 


3.24 


TO 




6050 


6056 


3.82 


3.82 


6.48 


6.49 


7.40 


7.22 


2.80 


1.08 


3.7 


3.7 


M67 


1.32 


6757 


6758 


4.36 


4.36 


2.95 


2.96 


11.88 


18.90 


3.18 


3.18 


4.26 10^^ 


4.26 10-^ 


ZAMS 




6476 


6499 


3.99 


4.00 


5.88 


5.79 


10.75 


9.05 


2.93 


1.23 


3.19 


3.11 


TO 




5824 


5724 


"ill 


3.75 


6.38 


6.20 


5.84 


5.73 


2.65 


0.93 


3.7 


3.7 


M67 


1.33 


6787 


6787 


4.36 


4.36 


3.06 


3.06 


12.09 


18.92 


3.17 


3.18 


3.88 10-- 


4.02 10-^ 


ZAMS 




6479 


6508 


3.97 


3.99 


6.11 


6.05 


10.71 


9.15 


2.92 


1.32 


3.14 


3.07 


TO 




5266 


5211 


3.70 


3.68 


4.98 


4.89 


3.63 


AAA 


1.82 


0.19 


3.7 


3.7 


M67 


1.34 


6814 


6808 


4.35 


4.35 


3.16 


3.17 


12.09 


18.97 


3.17 


3.17 


3.66 10-^ 


3.78 10-^ 


ZAMS 




6502 


6526 


3.97 


3.98 


6.28 


6.22 


10.71 


9.24 


2.93 


1.37 


3.05 


2.99 


TO 




5010 


4978 


3.54 


3.49 


5.88 


6.59 


2.54 


3.08 


1.44 


-0.15 


3.7 


3.7 


M67 


1.35 


6841 


6842 


4.35 


4.35 


3.28 


3.28 


12.09 


18.97 


3.17 


3.18 


3.38 lO-'' 


3.56 10-^ 


ZAMS 




6552 


6544 


3.98 


3.98 


6.43 


6.41 


10.90 


9.33 


2.93 


1.38 


2.94 


2.91 


TO 




4909 


4940 


3.29 


3.37 


9.86 


8.35 


1.88 


2.73 


1.30 


-0.23 


3.7 


3.7 


M67 


1.36 


6875 


6877 


4.34 


4.34 


3.34 


3.34 


12.23 


19.33 


3.17 


3.17 


2.95 10-- 


3.02 10-^ 


ZAMS 




6559 


6521 


3.97 


3.95 


6.63 


6.74 


10.81 


9.33 


2.94 


1.37 


2.84 


2.91 


TO 




4709 


4752 


2.78 


3.02 


26.90 


22.87 


1.11 


1.60 


1.26 


-0.37 


3.7 


3.7 


M67 


1.37 


6898 


6902 


4.34 


4.34 


3.51 


3.51 


12.23 


20.45 


3.17 


3.17 


3.03 10-^ 


2.98 10-^ 


ZAMS 




6515 


6555 


3.94 


3.95 


6.95 


6.88 


10.40 


9.87 


2.93 


1.15 


2.88 


2.82 


TO 




4637 


4453 


2.61 


2.23 


37.87 


76.04 


0.94 


0.91 


1.03 


-0.61 


3.7 


3.7 


M67 



here that were obtained by multiplying the values of Table |5] 
by I following Chandrasekhar (I1950I I. They mimic a mean 
V sin i associated to what can be considered a mean equatorial 
velocity given by stellar evolution models. The initial velocity 
and braking parameters (when relevant) were chosen to obtain 
surface velocities between 8 and 1 km.s"' at 3.7 Gyr For 
models with vzams — 19 km.s"' , the evolution of the surface 
rotation is driven by the radius changes as the star evolves, so 
that it essentially drops from the ZAMS to the upper RGB. The 
more massive of our models, which are also the more evolved, 
present the lowest surface rotation velocities, in agreement with 
the observations. The angular momentum lost by the wind at the 
evolutionary phases considered is accounted for, but is found to 
be negligible. 

Concerning the models that undergo magnetic braking, the 
figure shows the result of our calibration of the braking param- 
eter K. The extraction of angular momentum when using the 
Kawaler law is maximum when the model arrives on the ZAMS, 
because it is proportional to the surface angular velocity. By the 
time the model reaches the turn-off, it is efficiently slowed down 
and the surface rotation velocity has dropped by more than a 
factor of 2 as can be seen in Table |5] The evolution of v^^ii on 
the subgiant and red giant branches is then dominated by the 
increase of the radius as for the other family of models. 

Referring back to Talon & Charbonnel ( 120031 ). we note that 
our models at the age of the Hyades have effective tempera- 



tures between 6700 K and 6950 K, and would thus lie on the 
hot side of the Li dip. The convective envelope on the main se- 
quence is thus shallow, and internal gravity waves are not ex- 
pected to efficiently transport angular momentum or to influence 
the transport of chemicals. On the contrary, based on Palacios et 
al. (2003), we expect rotational mixing (e.g. meridional circula- 
tion and shear induced mixing) to efficiently transport Li so that 
the abundances derived from observations are reproduced. Let us 
note that the stars of M67 that are presently on the main sequence 
are less massive than 1.30 M©, and probably lie on the cool side 
of the Li dip. Thus according to Talon & Charbonnel ( 120031 ). we 
should expect an efficient transport of both angular momentum 
and chemicals by internal gravity waves, which should therefore 
be taken into account if these stars were to be modelled. 

3.3. Lithium as tracer of rotational mixing on ttie main 
sequence 

After chosing the masses and calibrating the initial rotation 
parameters so that the gravity and the vsini patterns as a 
function of effective temperature are well reproduced for our 
sample stars, we look at what is obtained for the lithium 
abundance. The results are presented in Fig. |7] where we plot 
the derived lithium abundances for our M67 sample and the 
abundances obtained for each of our models at 3.7 Gyr Li 
abundances appear here as a clear separator between the two 
families of models we considered. Models that are slow rotators 
and do not undergo magnetic braking on the ZAMS are ruled 
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Fig. 6. i/sini versus TgfF for our sample of stars. Symbols are 
the same as in Fig.|5] The theoretical predictions here are mean 
vsini and correspond to the Vgq values at 3.7 Gyr given in 
Table IS multiplied by |. The filled dot is for S1242, which has 
been identified as a peculiar Li-rich subgiant in a previous work 
of ours (Canto Martins et al. l2006l l. 
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Fig. 7. ALi versus Tefj for our sample of stars. Symbols are the 
same as in Figs.|5]and|6] 



out. Actually the Li abundance evolves similarly to standard 
models in these models and lithium barely decreases at the 
surface during the main-sequence evolution, as can be seen 
from Table |5] On the contrary, models arriving on the ZAMS 
with mild rotation (about 19 km.s"' ) and undergoing magnetic 



braking experience a reduction of the Li abundance by more 
than 1 .6 dex during the main-sequence evolution. Afterwards, 
when the models evolve to the RGB, the 1 " DUP further dilutes 
the Li of the envelope, and it almost disappears at the surface 
as can be seen from the 1.35, 1.36 and 1.37 M© models. These 
models, presented as filled pentagons on Fig. |7] follow the 
trend obtained for our M67 sample and nicely reproduce the 
ALi values derived from observations. 

Without braking the angular velocity gradient below the 
convective envelope is shallow when the star is on the main se- 
quence. As a consequence, the shear instability does not develop. 
Figure [8] presents the inner profile of Au and of the diffusion 
coefficient for chemicals (D - D„„r.circ. + Dshear + Dmicro) as a 
function of the radial coordinate in the 1.33 model while on 
main sequence (t 1.4 Gyr). The solid lines refer to the model 
without braking, and we see that the diffiision coeflicient is small 
in the region where Li is burned and the transport is strongly 
ineflicient. Because the meridional circulation velocity also 
remains very slow (a few 10"^ cm.s'' at most), the transport of 
chemicals is dominated by the microscopic diffusion, which is a 
slow process (D^ic ~ 100 cm^.s"') leading to a small decrease 
of the surface lithium abundance, which would also be obtained 
in a non-rotating model with microscopic diffusion. On the 
other hand, the extraction of angular momentum associated with 
the magnetic braking leads to an increase of the angular velocity 
gradient in the stellar interior in a way that the shear instability 
sets in and the meridional circulation becomes stronger as the 
star is slowed down on the ZAMS. For the 1.33 M© model, as 
can be seen in Fig. [8] the diffusion coefficient is two orders 
of magnitude higher than in the previous case just below the 
convective envelope, where the meridional circulation velocity 
is the highest (^ 10"^ cm.s"'). It is also more than one order of 
magnitude larger in the region where Li is destroyed. As a result, 
as can be seen from the Au profile, lithium (and hydrogen) 
diffuses inwards, reducing the gradient around the radial coordi- 
nate 0.65 Rq, and its surface abundance decreases. The depth of 
the DUP is not affected by rotation, and the subsequent decrease 
of ALi at this stage is similar, whether angular momentum was 
extracted or not during the early evolution on the main sequence. 

Let us emphasize that the lithium abundances obtained 
in our models with the adopted initial velocities and braking 
parameters fit the upper envelope of the values derived for the 
stars in our sample. The abundance determination lead only 
to upper limits for stars undergoing the DUP (in the range 
Tetf e [4700 K, 5600 K]). However, we do not expect the actual 
values to be much lower than these limits considering that for 
the stars with Teff ~ 4500K in our sample, Ali has been 
measured and is of about -0.5 dex. This value gives the lower 
limit for the post-DUP Ali value. Because the shape of the 
ALi vs Teff pattern is well reproduced by our models, a slight 
adjustment of the rotation parameters so as to decrease more the 
Li on the main sequence would concomitantly bring down the 
values for more evolved models. 

As explained above, we calibrated the braking parameters 
in our models to reproduce the evolution of the surface rotation 
velocity of stars beyond the turn-off in M67. This calibration de- 
pends on the equatorial rotation velocity adopted on the ZAMS. 
The models presented in Table |5] show that for VeqZAMS - 19 
km.s"' , the amount of angular momentum extraction needed 
to reach ysin/ ^ 8.5 km.s ' at the turn-off following Melo et 
al. (120011 ) induces an internal transport of chemicals that results 
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in a decrease of the surface lithium abundance that agrees very 
well with our observational data. Considering the uncertainties 
on the equatorial rotation velocities of the now giant stars of 
M67 when they arrived on the ZAMS, one may now wonder 
what would happen if we had chosen a faster rotation on the 
ZAMS. We thus computed models of 1.30, 1.33 and 1.36 Mq 
with a larger initial content of angular momentum leading to 
Veq,zAMS - 27.5 km.s"' , and calibrated the magnetic braking 
parameter of Kawaler's law in order to reach v sin / values in the 
observed range when the models reach 3.7 Gyr While the basic 
stellar parameters and the equatorial surface rotation velocities 
of these faster models resemble those presented in Table |5] the 
surface lithium abundance behaves very differently and is much 
lower at 3.7 Gyr This is the direct consequence of a stronger 
shear mixing that develops in the stellar interior because of the 
stronger magnetic torque applied at the surface of the star dur- 
ing its early evolution on the main sequence. This test allows us 
to constrain the ZAMS equatorial velocity of stars in the mass 
range 1.30 M© to 1.37 Mq of M67 considering that meridional 
circulation and shear turbulence are the main processes trans- 
porting lithium in this stars. We therefore exclude models with 
ZAMS equatorial velocities higher than 20 km.s ' . 

M = 1.33 Mg, X, = 0.45, [Fe/H] = -0.03 




Fig. 8. Internal profiles of Au (grey) and of the logarithm of the 
total diffusion coefficient (black) as a function of the radial co- 
ordinate in rotating models with initial mass of 1.33 M© with 
(dotted lines) and without (solid lines) braking. The profiles are 
shown for t ^ 1 .4 Gyr, when the models are on the main sequence 
(Xc - 0.46). Left ordinates refer to the scale for the diffusion co- 
efficient, while right ordinates represent the scale for the lithium 
abundance Au. Hatched regions represent convective zones. 



4. Conclusions 

We performed a spectroscopic analysis for a unique sample of 
27 evolved stars of the old open cluster M67 following an evo- 
lutionary sequence from the turn-off to the Red Giant Branch. 
In this analysis we used a spectral synthesis procedure based 



on MARCS model atmospheres to derive the stellar parameters 
(T(,yy, \ogg, {FejH}), microturbulent and rotational velocities, 
and the lithium abundance. The surface lithium content of the 
analysed stars follows a clear evolutionary pattern ranging from 
the turn-off to the Red Giant Branch. Our results reproduce 
the lithium abundance dispersion previously established for 
the M67 main-sequence stars and also confirm the well known 
gradual decrease of Au for Teff below 5500 K. The lithium 
abundances derived for our sample stars clearly indicate that 
this element has been already depleted at the stellar surface 
during the main-sequence stage, which also points to a depletion 
in excess of standard predictions. The present study largely 
extends the number of evolved stars studied in M67. For the 
first time, it offers solid possibilities for an evolutionary study 
of different stellar parameters {Au and v sin /) and also to test 
non-standard transport processes on homogeneous sample of 
turn-off and evolved stars of M67. 

We tested the hypothesis that lithium and surface rotation 
evolution might be connected by comparing our data to the pre- 
dictions of a set of stellar evolution models including the trans- 
port of matter and angular momentum by meridional circulation 
and shear-induced turbulence. It appears that these models (fol- 
lowing Maeder & Zahn 1 1 9981 formalism for the transport of an- 
gular momentum) consistently reproduce the entire set of ob- 
servational data when the evolution of angular momentum in 
the models includes an early phase of magnetic braking "a la 
Kawaler", the efficiency of which is calibrated on the v sin / data 
available for main-sequence and turn-off stars of M67. This re- 
sult is similar to what was also obtained by Smiljanic et al. 
(120101) for the open cluster IC 4651, and confirms the impor- 
tance of rotation-induced mixing in determining light elements 
abundance patterns in low-mass stars in particular for those with 
initial masses larger than 1.30 M© when at solar metallicity. 
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Table 6. Iron line list used for curve of growth and spectral s 
thesis, with atomic parameters, corrected oscillator strengths 
solar and Arcturus equivalent width measurements. 
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